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Abstract. We report the detection of the 2 P 3/2 -» 2 P m fine- structure line of C + at 157.74 yian in SDSSJ1 14816.64+525150.3 
(hereafter Jl 148+5251), the most distant known quasar, at z = 6.42, using the IRAM 30-meter telescope. This is the first 
detection of the [CII] line at high redshift, and also the first detection in a Hyperluminous Infrared Galaxy (L FIR > 10 13 L G ). 
The [CII] line is detected at a significance level of 8cr and has a luminosity of 4.4x 10 9 L Q . The L[ciij/Lfir ratio is 2x 10~ 4 , about 
an order of magnitude smaller than observed in local normal galaxies and similar to the ratio observed in local Ultraluminous 
Infrared Galaxies. The [CII] line luminosity indicates that the host galaxy of this quasar is undergoing an intense burst of 
star formation with an estimated rate of ss 3000M Q yr~'. The detection of C + in SDSS Jl 148+5251 suggests a significant 
enrichment of metals at z ~ 6 (age of the universe ~870 Myr), although the data are consistent with a reduced carbon to oxygen 
ratio as expected from chemical evolutionary models of the early phases of galaxy formation. 

Key words. Galaxies: high redshift - Galaxies: ISM - quasars: individual: Jl 14816. 64+525150.3 - Infrared: galaxies - 
Submillimeter - ISM: abundances 

1. Introduction In local galaxies, with far-infrared luminosities 1 Lfir < 

1O 12 L , the fractional contribution of the C + luminosity 
The 2 P V2 -> 2 P 1/2 fine-structure line of C + at 157.74 /jm t0 the far-infrared luminosity, which is a measure of the 
is often the brightest emission line in the spectrum of galax- E as heating efficiency, is typically in the range -3 < 
ies, ac counting for as much as -0.1-1% of their total lumi- log(L [C n] /L FIR ) < -2 (e.g., Stacey et al. 1991). For sources 
nosity dCrawford et allll985UStacev et allll99lHWright et all with far-infrared luminosities in excess of 10 12 L or with in- 



creasing star formation activity, this ratio drops b 


y an order of 


magnitude ( Malhotra et al., 


2001; Luhmanetal., 


1998, 2003; 


Negishi et all 


2001). Various possibilities have been proposed 



1991). The line is emitted predominantly by gas exposed to 
ultraviolet radiation in photo dissociation regions (PDRs) as- 
sociated with star forming activity (even in galaxies hosting 

AGNs), and has been extensively used to investigate the phys- t0 explain this effect at high luminosities, specifically: 1) a 

ical conditions of PDRs and to trace star formation in exter- high ra t io ° f ultraviolet flux to gas density, resulting in posi- 

nal galaxies (e . g. iMalhotra etafl 1200 it iLuhmati et all l2003| lively charged dust grains which in turn reduce the efficiency of 

Kaufman et allll999UGenzel & CesarskvN2000UBoselli eTall the gas heating by the photoelectric effect (e.g.|Kaufman etalj 

2002HPierini et alll2003l) . 1 1 99911 : 2) opacity effects which we aken the [CII1 emission 

line in infrared luminous galaxies fcferin A Phillinsl |2000; 
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1 We adopt the Lfir = L(40 - 500/jm) definition given in 
ISanders & Mirabell ll996l) . 
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iLuhman et ail ll998): 3) contribution to Lfir from dust asso- 
ciated with HII regions ( Luhman et al l 120031): 4 ) contribution 
to LpjR from an AGN (see Malh otraet all l2001). 

The possibility of probing the interstellar medium and trac- 
ing star formation in galaxies at cosmological distances by us- 
ing the bright C + emission line, when redshifte d into submil- 
limete r atmospheric windows, was propose d by Petrosian et alJ 
ill 9691). and further discussed bv lLoeblll99"l . lStarkllll997l) and 
Blai n et all J2000t) . However, [CII] has so far only been de- 
tected in local galaxies (z < 0.1) fro m space and airborne ob- 
servatories (e.g. lLuhman et alll2003l) . while all the searches for 
C + at higher redshifts have been unsuccessful. Deep searches 
for C + were carried out in z > 3 infrared luminous galaxies 
and quasars having massive reservoirs of molecular and neu- 
tral gas detected in CO emission lines. Specifically, C + was 
searched for in the z = 3. 13 damped Lyq ab sorption system 
towards PC 1643+4631 A Jlvison et all Il998l) . the z = 4.12 
quasar PSS 2322+1944 (Benford e t al., in prep.), the z = 4.69 
quasar BR 1202-0725 Jlsaak et all fl994t Ivan der Werd 1 19991 
Benford et al. in prep.), a z — 4.92 l ensed galaxy in the cluster 
CL 1352+62 jMarsden et all l2005l). and the z = 6.42 quasar 
SDSS Jl 148+5251 jBolatto et all l2004 the same object dis- 
cussed here). 

In this letter, we report the detection of C + in 
J114816.64+525150.3 (hereafter J1148+5251) the most dis- 
tant quasar known whose redshi ft of 6.42 puts it at the end 
of the epoch of re-ionization jFanetall EHoi. Millimetric 
and submillimetric continuum detections indicate that the host 
galaxy of this quasar radiates a far-IR luminosity 2 of * 2 x 
10 13 L o , which, if ascribed to starburst activity (and not to 
the QSO itself), would imply star formation rates of a few 
times 10 3 M yr 1 JBertoldi et alll2003atlRobson et alll2004l) . 
A massive reservoir of molecular and neutral gas of a few 
10 10 M(7i has been detected via CO line emission llWalter et all 
120031 120041 iBertoldi et all l2003bl) . The detection of copious 
amounts of dust and molecular gas indicates that the insterstel- 
lar medium in Jl 148+5251 is significantly enriched in metals. 
This enrichment is also seen through the detection of promi- 
nent ir on fine emission in th e rest-frame UV spectrum of the 
quasar jMaiolino et all \200A although the latter samples its 
nuclear region. 

The deep observations described in thi s paper, improv e sig- 
nificantly the noise-level reported by Bol atto et all ll2004) and 
enabled us to achieve the first detection of the C + emission 
line at cosmological distances as well as the first detection 
of this line in a galaxy with an infrared luminosity in excess 
of 1O 13 L . We discuss the implications of this detection for 
the physical properties of the ISM and of star formation in 
J1148+5251 3 . 



6.4 



Redshift 



6.12 



6.44 



2 Thanks to new submm data Beelen et al. (in prep.) obtained a new 
estimate of the far-IR luminosity, which is in the range 1.2 - 3.2 x 
10 13 L Q . 

3 In this letter we assume the concordance A-cosmology with Ho = 
71 km s-'Mpc- 1 , £1 A = 0.73 and fi m = 0.27 JSpergel et alll2003l) 
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Fig. 1. Spectrum of the [CII] 157.74 /urn emission line in the 
quasar Jl 148+5251 at z = 6.42 shown with a velocity resolu- 
tion of 56kms _1 (top panel) compared to the CO(6-5) emis- 
sion line (bottom panel - from Bertoldi et al. 2003b). The 
dashed curves show the gaussian fits to the line profiles (see 
Table[TJ. 



2. Observations and results 

Observations were carried out with the IRAM 30m telescope in 
two observing runs: February 2-3 and March 8-10, 2005. We 
used alternatively the C and D 150 receivers with filter banks 
covering a bandwidth of 1 GHz with 256 channels spaced by 
4 MHz. The receivers were tuned to 256. 1753 GHz correspond- 
ing to the frequency of the [CII] emission line (rest frequency 
1900.539 GHz) redshifted to z = 6.4189, the redshift of the 
CO emission lines, but we also observed with the frequency 
offset by ±100 km/s to minimize, and check, the effects of pos- 
sible instrumental artifacts. At this frequency the Half Power 
Beam Width is 9.6" and the 1GHz bandwidth corresponds to 
1170kms _1 . The observing conditions were generally good: 
low opacity, with T256GHZ < 0.1 during the March run and 
T256GHz < 0.2 during the February run, and system temperature 
T sys ~350 K. The observations were done in wobbler switch- 
ing mode with a switching frequency of 0.5 Hz and a wobbler 
throw of 200". In total, the source was observed for 12.4 hours. 
Calibration was obtained every 15 minutes using the standard 
hot/cold-load absorber measurements. Pointing was checked 
about every 1 .5 hour by means of one of the two 3mm receivers 
(A100/B100) and the pointing accuracy was found to be in the 
range 1 - 3". The focus was checked about every 3-4- hours, 
and was found not to change significantly. 

The data were processed using the CLASS software. After 
dropping some bad scans, only linear baselines were subtracted 
from individual spectra. The individual scans were realigned in 
frequency and co-added by weighting with 1/rms 2 . The result- 
ing profiles were re-gridded to a velocity resolution of 56 km/s 
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Tabl e 1. Properties o f the [CII] line observed toward SDSS Jl 148+5251 compared with the CO(6-5) line observed by 
Bert oldi et aiT ( l2003bl) . 



Line 


Vrest V obs 




Peak int. 


AV FWH M 


I 


L 




[GHz] 




[mJy] 


[kms^ 1 ] 


[Jy km s~'] 


[10 9 L G ] 


[cii] ep 3/2 - 1 Pi/ 2 ) 


1900.54 256.172 


6.4189+0.0006 


11.8 


350+50 


4.1±0.5 


4.4±0.5 


CO (6-5) 


692.473 93.204 


6.4189±0.0006 


2.45 


279+65 


0.73+0.076 


0.29+0.02 



leading to an rms of about 0.3 mK (2.8 mJy) 4 . We estimate the 
flux density scale to be accurate to about 30%. 

The final spectrum of the [CII] emission line in 
Jl 148+5251 is shown in Fig.^ The C + fine-structure line is 
detected with a confidence 8 cr. The velocity integrated flux is 
0.44 + 0.05 Kkm s _1 (in antenna temperature, T m b = 1.91 Ta), 
or 4.1 ± 0.5 Jykms -1 , and the line center corresponds within 
the uncertainties to the center of the CO line emission (Tabled 
and Fig. Q). 

Bola tto et alJ J2004I) used the James Clerk Maxwell 
Telescope to detect the C + line in Jl 148+5251, but could only 
place an upper limit to its line strength. To compare their 
limits with our result we note that the peak intensity of our 
line (12 mJy) is three times lower than the rms (32 mJy) in 
their final spectrum smoothed to a resolution of 10 MHz (or 
12kms _1 ). The [CII] velocity integrated flux reported in this 
paper is twice larger than the lcr limit estimated by Bolatto et 
al. for the intensity of any putative line in the JCMT spectrum. 



3. Discussion 

The emission of C + in Jl 148+5251 is an extreme example of 
what is seen in local infrared galaxies. The Ljchj/Lfir ratio in 
Jl 148+5251 is 2 x 10~ 4 , about an order of magnitude smaller 
than in local starburst galaxies. This is illustrated in Fig. |2] 
where we plot the Ljchj/Lfir ratio in nearby star-forming 
galaxies as a function of their far-infrared luminosity, together 
with the current upper limits for high-redshift sources and the 
value obtained for Jl 148+5251. A decrease of the L[ch]/Lfir 
ratio with increasing Lfir beyond 10 • L is clearly apparent 
in Fig. 13 The L[ch]/Lfir ratio for Jl 148+5251 is consistent 
with this general trend. Note also that our sensitive detection is 
consistent with the upper limits inferred for other high-z QSOs 
in Fig. |3 

The detection of strong [CII] emission suggests that the in- 
terstellar medium was already significantly enriched with met- 
als at z — 6.4, i.e. when the universe was as young as 870 Myr. 
However, the chemical evolutionary models predict that the en- 
richment of carbon relative to the a elements is delayed (since 
carbon is mostly produced by intermediate/low mass stars), and 
in particular the C/O ratio should be about 1/6 in the young 
stages of massive elliptical galaxies, while the absolute abun- 
dance of O should be simil ar to the local value, at least at an 
age of about 1 Gyr (Pipino & Matteucci, 2004). To achieve 
a better physical understanding of this source, and in partic- 
ular to investigate whether the C-poor scenario is consistent 



-2.5 



-3 



" -3.5 



° o 2 

o o°o°| ® 

g o 

o | o' 

o ° o 



o • 

M82 O 



x * 

I 



o ap.[CII] < lOkpc 
!<0.1 • ap.[CII] > lOkpc 
X ULIRGs 

4<z<5 ■ 

z=6.42 ■ 



10 



12 



14 



log L FIR (L ) 



4 The conversion between antenna temperature and flux density 
scale at 256 GHz for the 30m is 9.5 Jy/K. 



Fig. 2. The L[ciij/Lfir ratio versus Lfir for normal and 
starburst local galaxies (circles), as measured by the ISO 
Long Wavelength Spectrometer (LWS), and for high-z sources 
(squares), as measured from the ground. Jl 148+5251 is shown 
with a large square. The small squares indicate the upper limits 
for the three other high-z sources (identified by their names) 
where [CII] was searched for. Empty circles indicate galax- 
ies for which the aperture to measure [CII] (i.e. ISO-LWS) is 
smaller than 10 kpc, and therefore it may sample a smaller re- 
gion relative to the IRAS beam used for the far-IR flux. Crosses 
indicate local ULIRGs (Lir > 10 12 L G ). For reference, well- 
known local galaxies (Arp220, M82) are also identified. 



with the data or not, we have compared all available observa- 
tional data for the host galaxy of the quasar (i.e. [CII], FIR 
and CO transitions) with PDR models. The relevant observa- 
tional data are summariz ed in the second colum n of Tab|5] We 
used the PDR model of Kaufm an et alJ {l999), modified for 
the conditions appropriate for Jl 148+5251. In particular, we 
lowered the carbon-to-oxygen abundance as discussed above. 
We did not include detailed models for the evolution of dust 
jMorgan & Edmun ds! 120031). nor the p ossible contribution of 
dust from SNe ( Maioli no et all 120041) . we simply scaled the 
the abundance of dust grains and PAHs proportionally to the 
carbon abundance. The model requires high densities and high 
UV radiation fields to reproduce the data. However, a PDR 
model which fits observations of typical star forming regions 
fails to repro duce the strength of the high-J CO transitions ob- 
served by Bertoldi et al. (2003b). It is possible to increase the 
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Table 2. PDR model results (n = 10 5 cm" 3 , G = 10 38 , £ = 

1.8 x lcr'S- 1 ). 



Ratio 




Observed 


Model" 
C/0=l/6 C/0=local 


L[ciij/Lfir 




2.0(±0.9) x 10~ 4 


1.4 x 10~ 4 


3.8 x 10~ 4 


LcO(6-5)/L[CII] 




6.6(±2.0) x 1CT 2 


5.7 x lfr 2 


6.7 x 10- 2 


Lc0(7-6)/Lc0(6- 


-5) 


1.02(±0.17) 


0.85 


1.18 


Lc0(3-2)/Lc0(6- 


-5) 


0.22(±0.03) 


0.27 


0.17 


Lco(i-o)/Lco<6- 


-51 


< 2.4 x 1(T 2 


8.3 x 10~ 3 


4.4 x 10~ 3 



Note: " The abundances for the models with "local" ISM fractions are 
speci fically 0/H=3.2 10~ 4 and C/H=1.4 10 4 iSavase & Sembachl 
119961) . while for the "reduced" C/O model C/H=5.6 10" 5 and O/H has 
been kept to the same value as "local". 



gas heating in the CO-emitting region by increasing the cosmic 
ray ionization rate. Such a high ionization rate can mimic the 
higher gas heating due to an increased X-ray flux associated 
with the powerful QSO hosted by this system. We find a good 
fit to the observations with a gas density n = 10 5 cirT 3 , a FUV 
field Go = 10 3 8 times stronger than the average value in the 
Galactic ISM, C/0= 1/6, and a cosmic ray ionization rate f per 
H nucleus of ~ 10~ 16 s _1 , i.e. about 10 times the Galactic value. 
The cosmic background radiation at z — 6.4 has a temperature 
of about 22 K, much lower than the gas temperature observed 
in this system, it may however slightly affect the lower CO tran- 
sitions, but the latter are less relevant for our model (especially 
CO(1-0) which only has an upper limit). The resulting model 
reproduces all the observed line ratios and line-to-FIR ratios 
within a factor of 30%. This result shows that the C-poor sce- 
nario is consistent with the data. On the other hand it is not a 
proof of a low C-to-O ratio. Indeed, even using the C/O ratio 
typical of the local ISM the results of the model are still in fair 
agreement with the data, as shown in the third column of Tab.|3 
Summarizing, PDR models with high densities and enhanced 
X-ray (cosmic ray) flux can reproduce the observed properties 
of Jl 148+5251, both with local ISM abundances and with a 
reduced C/O ratio predicted by chemical evolutionary models. 
Tackling the issue of the chemical abundances in detail will 
require the observation of additional lines. 

As discussed in the introduction, the [CII] line is pro- 
duced by the UV field in star forming regions and therefore 
can be used as an indicator of the star formation rate (SFR). 
Bose lli et alJ J2002I) calibrate the luminosity of [CII] in terms 
of SFR, but their calibration only applies to low luminosity 
systems (Lfir < 10 10 5 L ). We have attempted to derive the 
L[cii] - SFR relation at high luminosities (Lfjr > 10 12 L o ). One 
possibility is to combine the relationship b etween far - IR lu- 
minosity and star formation rate derived by Kennicutt ( 1998) 
with the L[cii] /Lfir ratio predicted by the model (where we use 
the intermediate value of the Lpnj /Lfir obtained from the two 
models in Tab.2.), yielding the relation: 

SFR(M yr 1 ) * 6.5 lO" 7 (1) 

From Eq. 1 and from the [CII] luminosity we infer a star forma- 
tion rate in the host of Jl 148+5251 of SFR ~ 3000 M yr" 1 
consistent with the SFR estimated from Lftr jBertoldi et all 



2003a). Therefore our result demonstrates that the host galaxy 
of this quasar is the site of extremely powerful massive star for- 
mation, and it is probably tracing the earliest stages of the black 
hole-spheroid co-evolution. 

The first detection of the fine-structure line of C + in a high 
redshifted object shows the potential use of this emission line 
to investigate star formation and the physics of the interstellar 
medium in the early universe. Further detections and imaging 
of C + in high-z objects will continued to be made by using ex- 
isting instruments and forthcoming facilities, in particular with 
the Atacama Large Millimeter Array (ALMA) for which this 
topic is one of the primary scientific goals. 
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